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ABSTRACT 

Starless molecular cores are natural laboratories for interstellar molecular chemistry research. 
The chemistry of ices in such objects was investigated with a three-phase (gas, surface, and 
mantle) model. We considered the center part of five starless cores, with their physical con¬ 
ditions derived from observations. The ice chemistry of oxygen, nitrogen, sulfur, and complex 
organic molecules (COMs) was analyzed. We found that an ice-depth dimension, measured, 
e.g., in monolayers, is essential for modeling of chemistry in interstellar ices. Particularly, the 
H 20 :C 0 :C 02 :N 2 :NH 3 ice abundance ratio regulates the production and destruction of minor 
species. It is suggested that photodesorption during core collapse period is responsible for high 
abundance of interstellar H 2 O 2 and O 2 H, and other species synthesized on the surface. The 
calculated abundances of COMs in ice were compared to observed gas-phase values. Smaller ac¬ 
tivation barriers for CO and H 2 CO hydrogenation may help explain the production of a number 
of COMs. The observed abundance of methyl formate HCOOCH 3 could be reproduced with 
a Ikyr, 20K temperature spike. Possible desorption mechanisms, relevant for COMs, are gas 
turbulence (ice exposure to interstellar photons) or a weak shock within the cloud core (grain 
collisions). To reproduce the observed COM abundances with the present OD model, 1-10% of 
ice mass needs to be sublimated. We estimate that the lifetime for starless cores likely does not 
exceed IMyr. Taurus cores are likely to be younger than their counterparts in most other clouds. 

Subject headings: astrochemistry — molecular processes — ISM: clouds — ISM: abundances — ISM: 
individual objects (CB 17, CB 26, CB 27, B 68) 

1. Introduction 

Low-mass interstellar molecular clouds, such 
as the Bok globules, are a model case for cloud 
structure and star-formation studies, thanks to 

2009f) . For astrochemistry, stable starless cores 
within these clouds are particularly interesting. 

They are sufficiently dense for a significant mass of 
ice to accumulate onto the grains, supposedly not 
disturbed by gravitational contraction motions. 

Thus, modeling of surface chemistry can be sim¬ 
plified to calculations that consider constant, yet 
realistic physical conditions. 

The aim of the present paper is to study the 
evolution of the chemical composition of ices in 
starless cores with physical conditions that are 


their relative simplicity and proximity ( Hartmann 


relevant to actual cores. The modeling of sur¬ 
face chemistry in starless cores has a long history. 
Important papers that have presented significant 


19721; 

Allen & Robinson (1977): 

Pickles & Williams 

1977a 

bll: Tielens & Hagen (Il982 

): D’Hendecourt et al. 


Hasegawa et al. ( I992[l . A model for small, diffuse 


clouds (not the starless cor e) with s urface chem¬ 
istry has been presented bv iTurne'il ( 1998ll . Such 


two-phase models with constant physical param- 
eters have also been used in recent pape rs (e.g., 
Du et al. 2012 : Vasvunin fc Herb.^l20I3lJ l. 


Starless cores with physical conditions derived 
from actual obs e rvations have also been modeled 
( Du et al. 2012 : Maret et al.l 12013 : Lippok et al 
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2013t lAwad et alJ 12014 1. Recent advances in 


the field of interstellar ice chemistry in com- 
bination with availab le observational data (e.g. 
Launhardt et al. 2010ll make the modeling of ices 


in quiescent cores an interesting proposal. These 
advances include ices described as consisting of 


ing ayalaible for reactions 

(Hasegawa & Herbst 


1993b 

: Bergin & Langei 19 

97 I: Guppen & Herbst 


2007: 

Garrod & Paulvl 2011 

: iTaouet et al. I 2 OI 2 I). 


photoprocessing of subsurface ices (Kalyans & Shmeld 

2 OIOI: GarrodI 2013: Chang & HerbstI 2014). and a 


realistic model of ice molecule formation. 

An important aspect in interstellar ice models is 
a proper description of surface reactions in a con¬ 
tracting molecular core. N otably, this includes the 
surface synthesis of CO t (Ruffle fc HerbstJ 2001bl: 


Garrod fc Pauly! l2011ri . iBrown et al.l (jl988ri and 
Aikawa et al.l (j200lll show that the contraction 


phase is particularly important for a correct rep¬ 
resentation of the ice formation epoch. Recent 
dense core models that focus on surface chem¬ 
istry often do not consider the core contrac¬ 
tion phase fP u et al. 20121: Kalyans fc Shmeld 


2013 : Vasyunin fc Herbst 120131^ Chang fc Herbs! 
2014 : Reboussin et al.l 2014 ). The dependence 
of ice composition on the core density and in¬ 
terstellar extinction has be en shown by a num¬ 
ber of earlier studies (e.g., Rawlings et al.l 1992 : 
Viti fc WilliamsI 19991 : Green et al. 12001 ) 

We studied four starless cores - CB 17 - SMMl, 
CB 26 - SMM, CB 27 - SMM, and B 68 - with 
chemical kinetics modeling. They were treated as 
clumps of gas and dust with an increased density 
towards the center, surrounded by a more diffuse 
cloud. The cores undergo an initial contraction 
phase, until they reach a state with constant (ob- 
seryed) physical conditions (section 12.11) . A de¬ 
tailed surface and subsurface ice chemistry model 
has been used (j2.2|l . whose application to starless 
cores is the main noyelty of this paper. 

2. Methods 

The model ‘Alchemic-Venta’ was employed 
i n this researc h, taken from preyious works 
( KalyansI 2015alla . hereafter Papers I and II, re- 
spectiyely). It is base d on the ‘ALCHEMIC’ code 
( Semenoy et al.ir20I0l) . with details explained be¬ 
low. Physical parameters in the model were ad¬ 


justed for a selection of non-starforming cores. 
Papers I and II inyestigate collapsing prestellar 
cores, while the present paper focuses on stable, 
non-starforming cores. 


2.1. Physical model 

2.1.1. Cloud core conditions 

We inyestigate ice chemistry of four starless 
cores with their central density nu = n(H) -|- 
2n(H2), cm“^, and interstellar extinction Ay, 
mag, deriyed from obseryations. Only the core 
center part was considered (i.e., OD model), where 
riH and Ay are at their highest. Each simula¬ 
tion consists of a collapse Phase 1, when the cloud 
density increases, and a quiescent Phase 2 with 
constant physical conditions. 

A single-point OD model was chosen because the 
ices in the center of each core most likely do exist, 
despite a number of uncertain parameters. These 
include the efficiency of non-thermal desorp¬ 
tion mechanisms and turb ulent motions i n star ¬ 
less cores, as ob se ryed by iRedman et al.l (2006 ): 


Leyshakoy et al. ( 2014h , and Steinacker et al 


(120141) . Turbulence was not considered in the 
present study. Also, ices in core center likely haye 
experienced the full range of physical conditions of 
the core since its formation. The simple OD model 
puts certain constraints for the interpretation of 
results. 

Cores with a Plummer-like density profile 
( Plummeu llOllc IWhitworth fc Ward-Thomoson 
200 ll ). embedded in a surrounding cloud with a 
fixed column density A^^ut, were modeled. The 


Table 2: Initial abundances of chemical species. 


Species 

rii/rtH 

Ha 

0.4995 

H 

l.OOE-03 

He 

9.00E-02 

C 

1.40E-04 

N 

7.50E-05 

0 

3.20E-04 

Na 

2.00E-08 

Mg 

2.55E-06 

Si 

1.95E-06 

P 

7.63E-08 

S 

1.50E-06 

Cl 

1.40E-08 

Fe 

7.40E-07 
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Table 1: The constant physical parameters of modeled starless cores during their stable Phase 2. 


Core 

nH,o, cm~3 

Observational data 
ro, AU ri, AU 77 

Nout, cm 

Ref. 

Calculated data 

Ay, mag T, K ti, kyr 

CB 17 S - SMMl 

2.3E+5 

9.5E+3 

3.0E+4 

4.9 

0 


13.8 

8.2 

1465 

CB 17 L - SMM 

1.3E+5 

l.lE+4 

3.0E+4 

5.0 

3.6E+20 

Lb 

9.1 

9.2 

1451 

CB 26 - SMM 

8.7E+4 

8.0E+3 

4.0E+4 

3.0 

2.2E+20 

L 

6.5 

11.1 

1437 

CB 27 - SMM 

1.3E+5 

1.3E+4 

6.6E+4 

6.0 

3.1E+20 

L 

9.5 

8.9 

1451 

B 68 - SMM 

4.0E+5 

7.0E+3 

2.7E+4 

5.0 

3.0E+20 

L 

17.6 

7.9 

1482 


Note.— Assumed physical parameters: nH, 0 “density at the center of the core; ro—radius of the central density plateau; 77- 
power-law slope for density; r'l-radius of the region (core) in consideration; Aout-column density of the surrounding cloud. 
Calculated parameter ti is the length of the cloud core collapse Phase 1. 


^Schmalzl et al.l 11201411 
^Lippok et al.l 11201311 



Fig. 1.— Density of hydrogen atoms, interstellar extinction, and temperature of the B 68 SMM cloud core 
- an example of the evolution of physical conditions during the contraction Phase 1. After the contraction 
phase, these parameters were allowed to remain unchanged for a total integration time of 3.5Myr. 
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total column density to the core center can be 
calculated with 


1V„ = 




[1 + (r/ro)^]''/^ 


dr + lVout,cm (1) 


where nn.o is the (maximum) density at the center 
of the core, rg is the radius of the central density 
plateau, rj is the power-law slope at large radii, 
and ri is the radius of the core. The chemical 
modeling was done only for the center of the core, 
where r = 0 and the density profile of the core 
is necessary only for the calculation of Njj at its 
center. Table [I] summarizes the parameters uh.oj 
ro, ri, and fVout- The se are unique for ea ch core, 
and were derived from Lippok et al. ( 2013h for the 
CB 17, CB 26, CB 2 7, and B 68 cores and from 
Schmalzl et al. ( 2014h for the CB 17 core. The two 
‘versions’ of CB 17 are hereaf ter denoted CB 17 S 
and CB 17 L , for data from ISchmalzl et al.l and 
Lippok et al.l . respectively. This means that there 
are a total five models (CB 17 S, CB 17 L, CB 26, 
CB 27, and B 68) of four cores considered. 

Table m data are relevant to cores in the sta¬ 
ble Phase 2, when the contraction in Phase 1 has 
ended and the physical parameters of the cores 
are unchanging. This phase has an assumed in¬ 
tegration time of t 2 >2Myr for each modeled 
core. Prior to the stable phase, a core forma¬ 
tion phase was considered. The formation of star¬ 
less cores have been described before as a free-fall 


e.g.. iLee et al.l 120041) or ‘delayed free-fall’ (e.g.. 


Taauet et al.l 12014) . A ‘static’, step-like contrac¬ 


tion was employed by these authors. 

In the present study, the formation of the 
core was described as a del ayed free-fall p rocess . 
The approach pre sent ed bv IBrown et al.l (1198811 , 
Neiad et al. ( 199Clll and Rawlings et al. ( 1992 ) was 
used. A Plummer-like core with an initial peak 
density of nn.ini = 3000cm“^ contracts to the fi¬ 
nal density nu.o in a core formation time ti (i.e, 
length of Phase 1, see T able [J) . The parame¬ 


ter B in equation (1) of Rawlings et al. ( 19921) . 
which modifies the rate of the contraction, was 
taken to be 2/3. This means that the collapse to 
the final densities, indicated in Table [l] occurs in 
1.4-1.5Myrs. Figure [T] shows an example of the 
evolution of physical conditions for a core (B 68) 
during the contraction phase, up to l.SMyr. 

A v a lue B = Q.l has been u sed by 
( 199nl ). R.awlings et al. ( 1992 ). and 


Nejad^e^^ 


Ruffle et al 


( I 999 I) for prestellar core models. We take B = 
2/3 as a somewhat lower value for starless cores. 
The exact value of B has only a limited effect on 
ice composition, the main subject of the present 
study. This is because ice accumulation was cal¬ 
ibrated by modifying the poorly-known efficiency 
of indirect reactive desorption to match the ob¬ 
served proportions of H 2 O, CO, and CO 2 (sec¬ 
tion [2Tl4]). 

For the core formation epoch, rg was kept pro¬ 
portional to ( Keto fc CasellJ2(il()l : Taauet et al 
20141 ). The Plummer-like density profile of the 


contracting core was retained, while the density 
increases from UH.ini to nu.g- The value of ri, 
which determines the physical size of the region in 
consideration, was assumed constant. This means 
that the mass of the core (i.e., the number of 
H atoms within a sphere of radius ri) increases 
during the formation phase. A^out was retained 
constant, too. 

The the total integration time ti +12 was taken 
to be 3.5Myr and the results were analyzed for an 
interval of 1.0-2.3Myr (see section [37^ . In con¬ 
text with ice formation, the full time consists of 
three periods: (1) cloud with nH,o ~ 3 x lO^cm”^ 
and no ice on the grains; (2) ice accumulation up 
to a maximum thickness in a rapidly contracting 
core and up to 200kyr into the quiescent phase; 
(3) slow processing of ice by cosmic-ray induced 
photons after the freeze-out has ceased. The start 
of period (2) is marked with the appearance of 
first subsurface ice species at densities 5 x 10^ to 
1 X 10^cm“^. During the last period (3) the gas 
and the surface are in an approximate adsorption- 
desorption equilibrium. The nominal ice thickness 
decreases during period (3). 

Ice photoprocessing leads to a declining amount 
of CO in the icy mantles. Because of this, qui¬ 
escent Phase 2 integration times longer than 
2Myr were not considered. For longer times, 
the H 20 :C 0 :C 02 ice abundance ratio is no 
longer supported by observations of interstel- 


lar ices (iGibb et al. 

2004: Whittet et al. 

iBoogert et al. 20 111 

fiberg et al.1 20 111), f 


2007 


Mantle 


growth and the abundances of major ice species 
are considered in section [3] 

Gas and dust were assumed to be in thermal 
equilibrium with temperature deriv ed from the in- 
terstellar extinction Ay, in line with lGarrod fc Paulv 
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(|201lf) . In other words, it was assumed that the 
cores are heated externally by the interstellar ra¬ 
diation field. Ay was calculated with: 


Ay = 


JVh 


1.60 X 1021’ 


( 2 ) 


Table [T] shows the calculated values for Ay and T 
for each of the modeled cores in the stable phase. 

Table m shows the chemical abundances used 
at the start for all simulations. They are based 
on th e abu ndances derived by iGarrod fc Herbst 
( 2006ll and IWakelam fc Herbst ( 2008h . The rate 
of hydrogen ionization by cosmic rays was taken to 
be 1.3 X 10“^^s“^. The flux of cosmic-ray induced 
photons and interstellar photons was taken to be 
4875 and 10®cm“2s“i, respectively, with a grain 
albedo of 0.5. 

The self- and mutual-shielding of H 2 , CO, 
and N 2 was included wi th the use of the tabu¬ 
lated shielding functions ( Lee et aDll996l: Li et al 


2 OI 3 II . Shielding was also attributed to molecules 
in ice. It is important to note that the surrounding 
molecular cloud was included in these calculations 
by assuming an external H 2 column density of 
fVout/2. This was not included in Papers I and 
II that used a physical model built on similar 
principles. Such more accurate approach here is 
of particular importance, because more plausible 
values for indirect reactive desorption efficiency 
can now be derived (section 12. 1. 41) . The synthe¬ 
sis of hydrogen on grain surface was considered 
by, both, Langmiur-Hinshelwood and Eley-Rideal 
mechanisms, although the latter is of little impor¬ 
tance. 

2.1.2. Dust grain properties 

The grains were assumed to consist of two com¬ 
ponents. The first is inert, solid nuclei with a ra¬ 
dius of a = O.I/rm and an abundance 1.32 x 
relative to H atoms. In dark and dense conditions 
{Ay > 1.9mag), the second component, interstel¬ 
lar ice, begins to accumulate onto grain surfaces 
in significant amounts. The increase of grain size 
is calculated self-consistently, assuming a thick¬ 
ness of 350pm for each ML, which is multiplied 
by the number of ice MLs to obtain ice thickness 
b. The number of adsorption sites per ML was 
retained constant at 1.5 x 10®. Changes in this 
number likely do not influence t he agreement of 
model results with observations (Acharvva et al 


2011 ; Tag net et al.ll2012 ), and were not considered 
in the model. Dust grain albedo was assumed to 
be 0.5. 

The ice layer on the grains was described with 
the sublayer approach, developed in lPaoer^ . This 
means that a fully formed ice (accumulation from 
gas has ceased) consists of a surface layer (1-2 
MLs) and three subsurface layers - ‘sublayers’ - of 
approximately equal thickness. A sublayer nomi¬ 
nally may contain up to several tens of MLs. Two 
or more sublayers in a model allow to consider the 
depth-dependent composition of the ice mantles. 
Three sublayers w ere consid ered in the present 
study, similarly to Paper ij . For reference, they 
were numbered with the first being the shallowest 
and the third being the sublayer adjacent to grain 
nucleus. The transition from surface to the first 
sublayer and, subsequently, to second and third 
sublayers has been harmonized with the accretion 
rate, in line with Papers I and 11. 

2.1.3. Ice physical description 

The molecules are bound to the surface or into 
the mantle with their adsorption {Ed) or absorp¬ 
tion energies {Eb), respectively. Ed is the ‘basic 
value’, given in the reaction network, from which 
all the other energies were calculated. The binding 
energies, used for calculating the diffusion rates, 
for surface {E^^s) and mantle {E^^m) species were 
assumed to be O.dOAp and O.dOEs, respectively 
( Garrod fc Paulvl 2011 1. In Papers I and II it was 
argued that Ei,^M has to be higher, or similar to 
Ed, otherwise the movement of molecules in ice 
bulk would be faster than their evaporation from 
surface, which fits the description of a liquid, not 
solid. In line with the earlier works, it was as¬ 
sumed that Eb = ZEd, i.e., that typically E^^m = 
1.2Ed. For CO and CO 2 molecules, the ratios 
Eb,s/ED = Eh m /Er were taken to be 0.31 an d 
0.39, respectively ( Karssemeiier fc Cuppenll2014ll . 

In dark cores, the grains may be sufficiently 
cold to have a significant proportion of their sur¬ 
face covered by molecular hydrogen. H 2 provides a 
much weaker binding than water ice. The adsorp¬ 
tion e nergy for H and H 2 wa s calculated accord¬ 
ing to iGarrod fc Paulvl ( 201 ill . For other species, 
Ed was not affected by the surface coverage of H 2 
because of the large difference between the hop¬ 
ping rate of H? and vibration frequency of heavy 


molecules ( Paper Ilh . 
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Accretion of neutral gaseous species onto the 
grains with a radius a +6 was calculated with stick¬ 
ing probabilities of 0.33 for li ght species (H, H 2 ) 
and 1.0 for all other species (IKalvans fc Shmeld 
2 OIOII . A total of six desorption mechanisms 
were considered. The rate coefficients for evapo- 
ratio n and cosmic-ray-induced w hole-grain heat¬ 
ing ( Hasegawa fc Herbst 1993ah were calculated 
for each molecule as a thermal desorption at the 
equilibrium temperature and 70K, respectively. 
Desorption by interstellar and cosmic-ray-induced 
photons was considered with uniform yields for all 
species, 0.003 and 0.002, respectively. These val¬ 
ues were chosen to be in agreement with recent 
photodesorption experim ents with interstellar ice 
analogs, as discussed in Paper ij . Photodesorp¬ 
tion from subsurface mantle species was permit¬ 
ted, if the number of overlying MLs does not ex¬ 
ceed two. The yield for subsurface photodesorp¬ 
tion was reduced by a factor of 1/3. This approach 
was developed to be in line wit h the results of 
Andersson fc van Dishoeck ( 2008ll . 

A reaction-specific reactive desorption was 
used, wi th the prrk parameter set to 0.03 
( Paper ij l. This is the ratio of the surface-molecule 
bond frequency to the freq uency at which ene rgy 
is lost to the grain surface (IGarrod et al.ll2007ll . 


2.1.4-- Indirect reactive desorption 

Desorption induced by the energy released from 
the H-l-H exothermic surface reaction probably 
determines the relative proportions of major ice 
components H 2 O, CO, and CO 2 . In order to 
reproduce observations of interstellar ices, a se¬ 
quence CO > CO 2 >> H 2 O must be observed 
for indire ct reactive desorption efficiency /H 2 fd 
((Paper II I. Without indirect reactive desorption, 
the calculated abundance of solid carbon species 
in collapsing cores can be overestimated, when 
co mpared to the ob s ervational resu l ts pro vided 
by Gibb et al. (12004 h Whittet et al. ( 2007 1. and 
Oberg et al.l (|2011r) . For example, realistic CO and 
CO 2 relative abundances can be achieved, if the ice 
contains a su bstantial amount of CH 4 . H 2 CO, and 
(or) C H 3 OH ( Carrod fc Paulv 201ll : fTaauet et af 


2014) . Of the latter three species only methanol 
has been occasionall y observed in interstellar ices 
( Boogert et ahllioill ). 

The parameter /nafd can be described as the 
number of desorbed molecules per H 2 formed on a 


grain. Exact value for the efficiency are unknown, 
and, given the averaging and uncertainties in as- 
trochemical numerical simulations, it is reasonable 
to adjust /Hafd so that it produces an approximate 
fit with observations of starless cores. Such an 
approach can be used until a more precise value 
of /H 2 fd is found in experiments or, for example, 
quantum chemistry calculations. 

In order to derive /H 2 fd fo r the pres ent study, 
we employ a similar method to lPaper ij . A prestel- 
lar core in fee-fall gravitational collapse was con¬ 
sidered in that study. For the present research, 
the B 68 core model (Table [T]) was used. It has 
the highest final density and final Ay among all 
the starless cores in consideration. This allows 
the comparison of observed/modeled CO:H 20 and 
C 02 :H 20 ice abun dance ratios for hve A y values 
from the dataset of lWhittet et al.l ( 2007 ). An ad¬ 
ditional observational constraint for the model is 
the abundance of H 2 O, CO, and CO 2 ices at their 
respective threshold extinctions Ath- 

The above means t hat, while the B 68 cloud _ 

core i s currently stable ( Bergin et al.ll200'5lRedman et al 
20061 ) . we assume that it formed in a delayed free- 


fall process, which ceased when the core reached 
its present density. This assumption has been at¬ 
tributed to all simulated cores in this study. While 
it might not be entirely correct, it is probably just 
as useful as other realistic assumptions, because 
of the poorly constrained value of /H 2 fd, which 
has to be calibrated, regardless of the exact core 
formation path. 

Table 13] summarizes the observational and best- 
fit B 68 model results. The desorption efficiency 
/H 2 fd was calculate d according to the empiric ap¬ 
proach of jPaperll : 


/H2fd = Qexp{-Ejj/{Qi), 


(3) 


where Q = 9.5 x 10 ® and Qi = 10^. The adopted 
be st-fit valu es of Q and Oi are different from those 


Paper II . It is possible that the va lues obtai ned 


here are more reliable than those of Paper ij be¬ 
cause of the more advanced and realistic physical 
model (section r 2 . 1 .ip . 

The observational values for CO:H 20 and 
C 02 :H 20 ice abund ance ratios were taken from 
Whittet et al.l ( 2007 ). While newer datasets are 


available (Boogert et al.ll2011: lOberg et al.ll2011 ). 


IWhittet et al.l provides the Ay for each back¬ 

ground star, which enables to tie the observations 
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Table 3: Comparison of observed and calculat ed (B 68 core) CO and CO 2 ice abundances, relative to water. 
Observational data from IWhittet et al.l (|2007ll . unless noted otherwise. 


Source ID 

Av^ 

Observations 

CO a 0 / CO2 0 / 

HqO 

B 68 model 

CO 0 / CO2 0 / 

t'° , kyr 

043728.2-1-261024 

6.3 ± 1.5 


25.0 

1.8 

13.8 

1417 

042324.64-250009 

10.0 ± 0.5 

20.2 

19.1 

7.7 

13.0 

1451 

043325.94-261534 

11.7 ± 0.5 

12.0 

16.7 

12.5 

13.9 

1459 

043926.94-255259 

15.3 ± 0.5 

27.3 

16.7 

30.2 

16.5 

1471 

042630.74-243637 

17.8 ± 1.5 

45.1 

18.3 

43.2 

18.2 

1476 

Threshold ice abundances 


H 2 O, ML 

CO 2 , ML 

CO, ML 


Ath H 2 O 

3.2 it 0.1 


3.5 



1332 

Ath CO 2 

4.3 ± 1.0 



1.5 


1375 

Ath CO“ 

6.8 ± 1.6 




0.6 

1424 

Ath CO'i 

8.1 




1.3 

1437 


“Abundance ratio for icy species 


^’Integration time, corresponding to the 

IWhittet et all 1200 lb 

' iPaper 


Av 


value 


to specific time-points of the sim ulations with the 
corresponding Ay (see iPaner III , for more on this 
approach). 

2.2. Chemical model 

2.3. J. Reaction network 

Several changes have been introduced into 
the gas-grain network, which was acquir e d from 
Laas. Garrod. Herbst. fc Widicus Weaver (20111. 


The fou r new gas-phase reactions of IVasvtmin fc Herbst 
(I2ni3bll were added. All these reactions involve 
the methoxy radical CH 3 O. Because these authors 
did not discern between the methoxy radical and 
its isomer, the hydroxymethyl radical CH 2 OH, 
four complementary reactions were added for the 
latter species (1-4 in Table H]). 

The results of Paper ij show that formic acid 
HCOOH, formaldehyde H 2 CO, and formamide 
NH 2 CHO can be overproduced in subsurface man- 
tle conditi o ns w ith the surface reaction set of 
Laas et al. (I 2 OII), w hich has been adopted from 


Garrod et al.l (l2008fl . In an attempt to reduce 


their abundance, grain surface reactions 5-8 in Ta¬ 
ble [4] have been added to the network. They tend 
to reduce the production efficiency of HCOOH, 
H 2 CO, and NH 2 CHO because alternative outcome 
is now possible for their reactants. 

Reactions 5-7 were assumed to be barrierless. 
Reaction 8 has an assumed activation barrier 


i?yi=2100K, estimated from reactions of H 2 with 
compounds that have a similar structure to HCO. 
The barrier of the latter reaction is so high that it 
has little effect on the results of the current cold 
core model. 

Reactions 9-11 were added to more adequately 
reproduce the chemistry of sulfur oxides, with ac¬ 
tivation barriers taken from similar gas-phase re¬ 
actions. Finally, reaction 12 was added to avoid 
most of chlorine atoms ending up in solid CIO, 
which is probably unrealistic, given the availability 
of atomic hydrogen, irradiation, and an aqueous 
medium on grain surfaces. The activation barri¬ 
ers for four additional surface reactions - 13-16 in 
Table [4]- were also changed, in line with available 
data. 


2.2.2. Surface and bulk ice chemistry 


The rate for reactions on grain surface was cal¬ 
culated according t o the reaction-di f fusion com¬ 
petition approach of [Carroi^iPa^ ( 2011 ). The 
mo dified rate equations method has been applied 
Isee lPaper rt . As explained in Section [2.1.31 part 
of surface reaction products desorb into the gas. 
Reaction barriers are overcome either thermally 
by quantum tu nneling, whichever is faster 


or 


( Garrod et al.l[20(M . 


The rate of binary reactions in bulk ice was cal¬ 
culated for reactants that are immobile for most 
of the time (Papers I and H). It was assumed that 
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Table 4: Changes in reactions respective to the original network bv iLaas et al. ( 2011 1. 


No. 

New gas phase reactions®' 

a 

P 

7 

Reference 

1 

CH 2 OH + H — ^ 

CH 3 + OH 

1.6E-10 

0 

0 

NIST 

2 

CH 2 OH + O —^ 

H 2 CO + OH 

1.5E-10 

0 

0 

NIST 

3 

CH 2 OH + CH 3 - 

H 2 CO + CH 4 

1.4E-10 

0 

0 

NIST 

4 

CH 2 OH + CH 3 - 

C 2 H 5 OH 

l.OE-15 

-3.00 

0 

Vasvunin & Herbst ('20I3bl*’ 

No. 

New surface (ice' 

) reactions 

Ea,K 



Reference 


5 H + HCO —^ CO + H2 0 

6 OH + HCO—^CO + HzO 0 

7 NH2 + HCO —5- CO + NHa 0 

8 H2 + HCO —^ CH3O 2100'= 

9 OH + SO—^S02 + H 0 

10 O + HS —^ S + OH 956 

11 0 + HS—^SO + H 0 

12 H + CL—^HCL 0 


This work 
This work 
This work 
This work 
Gas phase 
Gas phase 
Gas phase 
Gas phase 


No. 

Changed surface (ice) reactions 

Ea, K 

Reference 

13 

CO -h 0 ^ CO 2 

290 

Roser et al. ('20011 

14 

02 -fH^ 02 H 

600 

Du et al. ('20121 

15 

H -h H 2 CO —^ CH 3 O 

2100 

'WooTi 120021 

16 

O 3 -f H ^ O 2 -f OH 

0 

Mokrane et al. 120091 


“In addition to the four reactions from I Vasvunin fc Herbs'^ ll2013ljl . Reaction rate coefficient k = a(T/300)^exp(— 7 /T) 

^Data adopted from CH 3 O + CH 3 reaction 

“^Estimate 




























each molecule vibrates in a lattice cell with ten 
neighboring molecules. Reactants have to over¬ 
come a certain energy barrier Eprox = O-IEd in 
order to achieve sufficient proximity for a reaction 
to occur. Molecule diffusion in ice and reaction- 
diffusion competition have been taken into account 
in a manner, similar to that for surface reactions. 
I refer the reader to Papers I and II for a more 
detailed justification and formalism for the calcu¬ 
lation of reaction rate coefficients. 


Surface and bulk-ice molecules are subjected to 
photodissociation by interstellar and cosmic-ray- 
induced photons. For surface species, the dissocia¬ 
tion rate has been taken to be equal to that of gas- 
phase species. For molecules within the mantle, 
the photon flux is attenuated by overlying MLs. 
Each ML has an assumed absorption pro bability 
of 0.007 ( Andersson fc van Dishoeckll2008l l. In or¬ 
der to calculate the photon flux, respective to each 
sublayer, the value rel evant to m iddle ML of that 
sublayer was used (see IPaoer II ). 


3. Results: ice formation and major com¬ 
ponents 

To avoid repetition in the studies of five roughly 
similar models of objects of the same class, the 
modeling results are not reviewed for each case 
separately. Instead, we focus on several astro- 
chemical problems, relevant for ices in starless 
cores. The use of five models with physical condi¬ 
tions derived from observations instead of an ab¬ 
stract model adds to the credibility of the research. 
We also have to keep in mind that the OD model 
adequately represents only the center part of each 
core. The abundances values used here are calcu¬ 
lated with N(A)/N(H), i.e., relative to total hy¬ 
drogen. This means that observationally detected 
abundances that are given relative to H 2 were di¬ 
vided by two for comparison with calculation re¬ 
sults. 


3.1. The evolution of the ice layer 

Ice layer is basically nonexistent (around 
O.IML) when Ay < 1.9mag, which means that 
the calculation results (abundances of species) are 
irrelevant for ice chemistry if t < IMyr. Because 
of this, we interpret results for integration times 
t > IMyr, only. Ice formation observes the rules 
outlined in IPaoer III . As the cloud contracts and 


Ay increases over 2 mag, ice molecules begin to 
accumulate on grain surfaces, forming a mixture 
of H 2 O, CO 2 , and NH 3 . The ice is dominated by 
interstellar photons, which means that all surface 
CO is oxidized to CO 2 via the sequence: 


H20(H0:0H(C0)C02, (4) 


probably the most important chemical transfor¬ 
mation in bulk ice. When the cloud becomes 
darker, the flux of interstellar photons is dimin¬ 
ished, rapid ice photoprocessing becomes impossi¬ 
ble, and CO begins to accumulate at Ay > 5. The 
accretion of CO is faster than the surface synthesis 
of CO 2 , and CO becomes the second most abun¬ 
dant ice molecule when Ay > 12mag. 

Figure [2] shows the nominal thickness of ice and 
the abundances in MLs for a few most important 
species in each sublayer for all five starless cores. 
Ice thickness peaks shortly after the end of the core 
collapse stage, and then slowly decreases. The rea¬ 
son for this decrease is the conversion of simple ice 
species into more complex ones via photoprocess¬ 
ing. Most importantly, the sequence O converts 
water and carbon monoxide into carbon dioxide. 
Thus, the number of molecules per grain and the 
nominal ice thickness are reduced. This can be 
seen in figure [3l A significant amount of hydrogen 
is released in the process. H 2 partially accumu¬ 
lates in the mantle when T < IQK. The inner 
sublayer 3 is the least affected because it has a 
lower proportion of CO and it is better shielded 
than the above l ayers. These results ar e differ¬ 
ent from those of Chang fc HerbstI ( 2014h thanks 
to the gradual increase of gas density and Ay 
during Phase 1, which ensure an ice composition 
that is much more differentiated between layers. 
Such a chemical differentiation (H 2 O in inner and 
CO in outer layers) has been confirmed by mod¬ 
els that consider ice forma tion in collapsing cores 


with monolayer ac curacy (|Garrod fc Paulyll2011 


Taauet et al.ll2014l) . 


CB 26 is the smallest of the fiye cores, it has 
the lowest final Ay and density yalues (Table [IJ. 
This means that molecules on grains deposit at 
a lower rate, and haye undergone significant sub¬ 
surface photoprocessing eyen before the freeze-out 
has ended. With the ‘freeze-out epoch’ we under¬ 
stand the interyal between the time of formation 
for the first sublayer (tsub) and the time, when 
maximum ice thickness is reached (tmaxi as speci- 
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1.0 1.5 2.0 2.5 3.0 3.5 


Fig. 2.— Calculated ice thickness and abundance in MLs for major species in the sublayers for the modeled 
cloud cores. The number in parenthesis indicates the number of the respective sublayer, while ‘Other’ stands 
for all other icy species in that sublayer. The vertical gray line indicates tso for each model. 
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e to H atom number density, of major species in ice for the five 
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Table 5: Calculated ice composition parameters in the starless core models. 


Core 

Ay 

^sub^ 

Myr 

ML 

Max. ice 
^max 5 Myr 

thickness 

CO 2 % 

tbO, 

Myr 

ML 

Ice at 3.5Myr 

CO cy CO 2 cy 

H,0’ H,0> 

CB 17 S 

13.8 

1.25 

170 

1.57 

81 

23 

2.31 

147 

68 

81 

CB 17 L 

9.1 

1.21 

164 

1.62 

74 

35 

1.92 

141 

61 

107 

CB 26 

6.5 

1.21 

153 

1.68 

69 

59 

1.57 

133 

48 

150 

CB 27 

9.5 

1.06 

158 

1.60 

69 

51 

1.58 

131 

51 

175 

B 68 

17.6 

1.15 

171 

1.52 

78 

25 

2.12 

146 

63 

90 


Note.— tsub~time for ice thickness to reach IML; tmax—time of maximum ice thickness; tso—time for the C02:H20 ice 
abundance ratio to reach 50 %. 


fied in tabled. ^50 is the point in time when the 
C 02 :H 20 ratio in ice has grown to 50%, and for 
CB 26 it occurs even before tmax- 

On the other hand, CB 27 is the largest core 
considered. Paradoxically, this results in a sim¬ 
ilar ice composition to that of CB 26. The Ay 
threshold for accumulation of a significant mass 
of ice on the grains (>1ML) is close to 1.9mag. 
This occurs at a density of ~ 5 x 10^cm“^ for 
CB 27, and in the range (8 — 9) x 10^cm“^ for 
the other four cores. The result is that the first 
ice sublayers in CB 27 appear some 150kyr earlier 
than in other cores. The formation of sublayer 1 
in CB 27 is complete some 20kyr earlier than for 
the other cores. This means that sublayer 1 has 
formed in a relatively diffuse environment, which 
is still largely dominated by interstellar photons, 
where water photodissociation products combine 
with CO on the surface to produce an ice that 
contains a large proportion of CO 2 . 

Because indirect reactive desorption, which 
largely governs ice composition, has been cali¬ 
brated for the B 68 core (section 12.1.41) , the ex¬ 
plained mechanisms may lead to different results, 
if /uafd is different. The important thing is that 
both, a diffuse medium or early ice formation, 
may lead to a significantly photoprocessed ice. 
The first is characteristic for low-mass molecular 
cloud cores, while the second is more likely for 
relatively massive cores. 

Table [S] summarizes data for each model at two 
important points in time: when ice thickness is 
at its maximum, and at the end of the simulation 
run. The thickest ice layer is attained in CB 17 S 
and B 68 models. High ice thickness for B 68 


can be expected because of its high extinction and 
final density. For CB 17 S, neither its density, nor 
extinction are the highest. A thick ice layer for this 
core arises because JVout for this model is 0 , which 
means that interstellar photons are attenuated to 
a lesser extent and the ice layer begins to form at a 
later, denser stage of the collapsing core. In turn, 
ice accumulates rapidly and undergoes less intense 
photoprocessing, and a smaller number of CO and 
H 2 O molecules are converted to CO 2 . This results 
in higher nominal ice thickness. This conclusion 
is supported by the abundance of carbon oxides 
- a higher proportion of CO 2 means a nominally 
thinner ice (Table [5]). 


The above analysis confirms that extinction 
largely regulates the formation of the ice layer 
on the grains because Ay determines the num- 
ber of interstellar photon s that reach the grains 
(|Watson fc SalDeteiill972ll . Interstellar extinction 
also affects the abundance of radical species on 
surfaces and in the gas, which determines the ef¬ 
ficiency of direct and indirect reactive desorption. 
It can be safely said that the formation history 
of a starless or prestellar core determines its ice 
composition. 


3.2. General ice chemistry 

Aside water and carbon oxides, the nitrogen 
species NH3 and N2 also can be considered major 
ice constituents. Their abundance is little affected 
by ice photoprocessing. Other ice species whose 
abundance may exceed 1% relative to water for 
at least one of the modeled cores include H2, O2, 
H2O2, O2H, and H2S. 

Figure [ 5 ] shows ice evolution up to 3 . 5 Myr, 
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the full length of the simulation runs. Although 
the modeled cores have different physical param¬ 
eters, ice sublayers share a similar composition in 
all cases. The inner sublayer 3 contains a large 
amount of water, in addition to CO 2 and NH 3 . All 
CO molecules here are quickly consumed, which 
makes chemical processes in sublayer 3 somewhat 
different to those in above layers. The lack of CO 
and its daughter atomic C means that O, OH, N, 
NH, S, and other radicals are more available. 

Sublayer 2 initially forms as a Ril:l mixture of 
H 2 O and CO, with additions of CO 2 and N 2 . As 
the time progresses, CO and H 2 O combine, and 
the sublayer becomes increasingly dominated by 
CO 2 . The outer sublayer 1 largely consists of 
CO, with an admixture of H 2 O, CO 2 , and N 2 
at roughly similar proportions. CO remains the 
main species in this layer, although the propor¬ 
tion of CO 2 is steadily growing. Sublayer compo¬ 
sition may profoundly affect the chemistry of mi¬ 
nor species because the main reacting species are 
photodissociation fragments of major molecules. 


3.3. Comparison with observations and 
age limits 


The initial chemical composition of ices is 
largely regulated by indirect reactive desorption, 
calibrated with the B 68 model (section 12.1.41) . 
Table m shows that carbon oxides are somewhat 
depleted during the earlier stages in cloud con- 
traction, compared to the observational data by 
Whittet et al.l ( 2007l l. We found it impossible to 


remove this discrepancy with adjustments in the 
efficiency of indirect reactive desorption. Equa¬ 
tion ®. This depletion may arise, for example, if 
the calculated abundance of atomic H is higher by 
a factor of two than H abundance in the real core. 
Another possible cause can be changes in sur¬ 
face properties as proportionally more and more 
CO molecules stick to the grains during the ice 
formation epoch. We will investigate the latter 
suggestion in a separate paper. 


Figure |4] shows that molecules, typically formed 
on grain surfaces or dense gas (e.g., H 2 O, CO 2 , 
NH 3 , and O 2 ), are also abundant in the gas phase 
during or immediately after the core contraction 
stage, when density is in excess of « 10 ^cm“^. 
Gas and surface molecules are actively inter¬ 
changed during the freeze-out epoch. Such gas- 
phase abundance peaks in a layer around a collaps¬ 


ing co re have been found in models of prestellar 
cores (Carrod fc Herbst 2006: Vasvunin fc Herbst 


2013al ; IPaoer Hr ) and observ ations of water vapo r 


in the prestellar core LI544 ( Caselli et al.l 2012fl . 
The present model predicts is that such peaks are 
likely to be shorter and higher thanks to more ef¬ 
ficient desorption mechanisms. Desorption retains 
a substantial amount of refractory species in the 
gas phase until very late times in core contraction. 

We can conclude that the Phase 1 modeling re¬ 
sults confirm that gas falling into a prestellar or 
protostellar core is rich with molecules synthesized 
on the surface. The abundance peaks predicted by 
the present study are higher than in other similar 
models, thanks to more efficient desorption mech¬ 
anisms. The main desorp tion agent are interstel lar 


photons, as suggested by Dominik et al. ( 2005ll . 


Subsurface ice processing may place constraints 
for the maximum age of starless cores. The ob¬ 
served solid C 02 :H 20 abundan ce ratio does not 
exceed 44% ( Booeert et al.l[20llh . although higher 
values are permitted by upper limits. As a thresh¬ 
old value for the age of the cores we chose the 
longest time taken for any of the cores to reach 
a C 02 :H 20 ice abundance ratio of 50%, corre¬ 
sponding to a time For the CB 17 S core 
tso = 2.31Myr, longer than t^o for other cores (Ta¬ 
ble [5|). A single end time of 2.31Myr was applied 
for all simulations. The decision to adopt a sin¬ 
gle and rather long total integration time for all 
models results in that for the CB 26 and CB 27 
cores C 02 :H 2 0 is higher than 50% for most of the 
time. With this approach the model considers a 
variety of icy environments and allows for an easy 
comparison between the five simulation results. 

The above means that the actual lifetime ^2 of 
a quiescent Phase 2 core with constant physical 
conditions is below O.OMyr for all models. This 
is significantly less than the initial assumption of 
t 2 R:!2Myr for a 3.5Myr total simulation length. 
Thus, the modeling results were analyzed for an 
integration time interval of 1.00-2.31Myr. 


If gas and dust in the center part of a starless 
core is not effectively mixed with its surroundings, 
C 02 -rich ices should be present in starless cores. 
Turbulence and mixing of gas in the dense cores 
is a poorly-known parameter t hat may contribute 


to the destruction of the ices (Boland & de .Tone 

1982 

: IWilliams Kr. HartauistI 1984: R.edman et al. 

2006 

: lAdams & Shul 20071: Levshakov et al.l 20141: 
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Steinacker et al. 2014ll . Given the uncertainties, 


we treat the chemistry of C02-rich ices as a theo¬ 
retical possibility. Their observational detection 
or non-detection might provide additional con¬ 
straints on the physico-chemical processes that 
govern ice formation, chemical processing, and 
desorption. 


4. Results: elemental chemistry 

4.1. Oxygen chemistry in ice 

4.1.1. Chemical processing of minor oxygen 
species 

Three major ice species - H 2 O, CO, and CO 2 
- contain copious amounts of the eighth element, 
and the chemistry of oxygen in ice is closely con¬ 
nected to these molecules. When water and car¬ 
bon oxides are excluded, the bulk ice chemistry 
of oxygen for CB 17 S and B 68 is initially dom¬ 
inated by molecular oxygen O 2 , which is rather 
quickly converted into hydrogen peroxide H 2 O 2 , 
as shown in figure [5l For other models, where 
ice photoprocessing is more intense, H 2 O 2 is the 
more important species from the beginning. H 2 O 2 
is more stable in irradiated interstellar ices than 
O 2 , which mainly arrives from the gas phase. The 
hydroperoxyl radical O 2 H is the most abundant 
radical for oxygen chemistry in ice, thanks to its 
high Eb,M. It is generated by two mechanisms - 
H addition to O 2 (early times) and the photodis¬ 
sociation of H 2 O 2 (late times). Hydroxyl radical 
OH begins to accumulate to a limited extent in 
the inner sublayer 3 only after the CO molecule 
has been oxidized. 

Interestingly, CB 26, which has the lowest Ay 
and highest flux of interstellar photons, has also 
the lowest radical content. This is because the 
higher temperature in CB 26 (see Table [T]) allows 
a greater radical mobility (especially for atomic 
H) in the mantle and, thence, reactivity. Low rad¬ 
ical content in CB 26 ices is characteristic also for 
other compound classes (sections 14.2114.3115.1|) . 

Because ice in CB 17 S is initially rich in O 2 (fig¬ 
ure [5]), it has also a higher abundance for H 2 O 2 , 
which forms via reaction 14 (Table H]). Substan¬ 
tial amounts of H 2 O 2 (relative abundance XH 2 O 2 
in excess of 10“®) in sublayer 3 for CB 17 S and 
B 68 are able to form from water photodissocia¬ 
tion products that combine together after all CO 


has been oxidized in that sublayer. 


4.1.2. Comparison with observations 

Figure [7] shows that the calculated abundance 
of gaseous hydrogen peroxide H 2 O 2 can be rea¬ 
sonably high for a short period, up to 9 x 10“^® 
relative to total hydrogen atomic number density 
for CB 17 S. The abundance of its associated 
O 2 H radical, however, is lower by three orders 
of magnitude. Both these species have been ob¬ 
served in interstellar mediu m with roughly simi- 
lar abundances of5xl0“^^ ( Bergman et al ] l2mi[ 


Parise et al.l 1201211 . It has been suggested that 


reactive desorption is responsible for tr ansferring 


these species from the surface to the gas (iDu et al 
20121 1. However, reactive desorption may be too 


inefficient to maintain a s ignificant amount of such 
heavy species in the gas ( Paper l3 l. The compari¬ 
son of abundances here is indicative, only, because 
of the limits of the OD model. 

Because practically all of H 2 O 2 in ice resides in 
the sublayers, photodesorption from subsurface ice 
may play a role in the existence of gas-phase H 2 O 2 
and O 2 H. More often than not, dissociative pho¬ 
todesorption occurs, when molec ule fragments are 
ejected, instead of intact species (lAndersson et al 


20061 : Andersson fc van Dishoeckll2008 1. This as¬ 
pect has not been c onsidered i n the present model, 
or in the model bv iDu et al.l (|2012ll . It likely that 
photodesorption of H 2 O 2 also induces a high gas- 
phase abundance of O 2 H. 

Beginning with Goldsmith et al.l (2002), molec¬ 
ular oxygen O 2 has been repeatedly observed in 
shocked gas. A few detections have been associ¬ 
ated with the cold dense core p Ophiuchi A, with 
inferred abundances in the range 2.5x 10“^...1.25x 
10“® relative to total hydrogen. Such a high 
abundances likely can be observed only during a 
short period in the core’s evolution ( Larsson et al.l 
2007; Liseau et al.l I2OIC1II . These results agree 
with the calculation results shown in figure [71 
where the maximum abundances lie in the range 
7 x 10“^...4.5 X 10“®. A quantitative comparison 
is possible in this particular case because of the 
centrally localized nature of the observed source 
( Liseau et al.l 2OIC1II . 
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Fig. 5.— Calculated abundance of selected oxygen speci( 
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Fig. 6.— Calculated relative abundance of main O and 0-H species in ice mantle sublayers for CB 17 S, 
and CB 26 models. The latter represents a significantly photoprocessed ice. The irregularities for CB 17 S 
surface layer at t « 1.5Myr are artifacts. 
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4.2. Nitrogen chemistry in ice 

4-2.1. Chemistry of important species 

Figures [5] and [5] show that ammonia NH 3 and 
molecular nitrogen N 2 are the main nitrogen reser¬ 
voirs in ice and gas. Nitrogen chemistry is regu¬ 
lated by proportions of N 2 and NH 3 in the surface 
and the three sublayers (section [321) • The rela¬ 
tively clear division of NH 3 (in sublayer 3) and N 2 
(in sublayers 1 and 2) means that N, NH, and NH 2 
radicals are not equally available in all sublayers. 
Interesting examples of relatively abundant minor 
nitrogen species in ice are formamide NH 2 CHO 
(discussed in section |S]), hydrogen cyanide HCN 
and isocyanide HNC. 

The chemistry of HCN and HNC is controlled 
by the proportions of CO, N 2 , and NH 3 in each 
sublayer. HNC is more favored by ice chemistry 
at earlier times than cyanide HCN, while HCN is 
the more favored isomer during late times. The 
two isomers form via different pathways. HNC is 
mainly synthesized in the reaction NH 2 -I-C, while 
HCN is the product of H addition to CN. The lat¬ 
ter forms via C-l-N. This means that HNC is a 
product of NH 3 and CO mixture photoprocessing, 
while HCN is synthesized in a N 2 and CO mix¬ 
ture. The HNC pathway is more efficient because 
it requires no subsequent hydrogenation. The re¬ 
sult is that 60 (for CB 17 S) to 90% (for CB 26) 
of both isomers are concentrated in the outer sub¬ 
layer 1, where CO and N 2 are highly abundant, 
and NH 3 is available, too. The remaining 10-40% 
reside in sublayer 2, while sublayer 3 is very poor 
on HCN and HNC (<2%), because of a lack of CO 
there. Because the surface closely interacts with 
gas-phase CO and N 2 , HCN is the main isomer for 
the top ice layer. 


4.2.2. Comparison with observations 


Recent observed NH 3 gas-phase abundances, 
lie in th e range 9 x 10~^. . .7.5 x 10~^ relative to 


total H (jRuosfeiienet^ 


20 lit Levshakov et al 


20141 1. ICrapsi et al.l (|2007l l derives a central NH 3 
abundance of 4.5 x 10“® for the L 1544 dense 
core. These values qualitatively agree with cal¬ 
culation results for cores with ages in the range 
t = 1.5...1.8Myr. 

Figure [10] shows the calculated HNC:HCN 
abundance ratio the center of the respective dark 


cores in gas and solid phases. We note that the 
observed gas-phase values of this ratio are in the 
range of 0.54-4.5, with no distinction between 


YSOs and starless cores ( Hirota et al. 1998ll . 


HNC:HCN ratio in ice may exceed unity for all 
cores, except for CB 26, where lower abundance 
of CO prevents the build-up of HNC. The highest 
calculated HNC:HCN ratio is in CB 27 ice phase, 
where it reaches 1.75 briefly at 1.54Myr. Such a 
peak arises because of effective HNC production 
in sublayer 1 , not observed in other simulations. 
In turn, this occurs because sublayer 1 in CB 27 
is richer in NH 3 by a factor of >2 than in other 
models. Finally, NH 3 is more abundant because 
ices in CB 27 are accumulated at lower extinc¬ 
tions, which means a higher proportion of species 
that form from free atoms on the surface (H 2 O, 
NH 3 , CH 4 , etc.). 


4.3. Sulfur chemistry in ice 

Figures[TT]and[T 2 ]show that H 2 S and SO are the 
most abundant sulfur species in ice and gas. H 2 S is 
more abundant than SO for most of the time. An 
exception is the CB 26 core, where H 2 S is the only 
dominant species and SO ice forms with a maxi¬ 
mum relative abundance of 5 x 10“® at 1.56Myr, 
only, thanks to surface reactions. Then, SO is 
quickly consumed by ice mantle reactions because 
of the higher temperature in the CB 26 core. In 
the other four cores, the HS radical accumulates 
and has time to react with atomic O, producing 
SO. Increased mobility of hydrogen in CB 26 ice 
does not allow the accumulation of HS, S, and 
many other radicals. H 2 S can be expected to 
be the main sulfur reservoir for photon-dominated 
ices. Notably, CB 26 still has a maximum SO 
gas-phase abundance comparable to that in other 
modeled cores (figure [T 2 |) . 

The abundances of minor sulfur molecules are 
affected by the proportions of major ice con¬ 
stituents. Molecules, whose synthesis is associ¬ 
ated with CO and its daughter radicals (primarily 
atomic C, atomic O is available also from CO 2 ), 
include OCS, CS, H 2 CS, and C 2 S. These are pri¬ 
marily associated with sublayers 1 and 2 (see fig¬ 
ure dsi). Oxidized forms SO 2 and NS arise in sub¬ 
layer 3. The above is true for cores CB 17 L, 
CB 17 S, CB 27, and B 68 . For CB 26, the di¬ 
versity is largely limited by the dominating H 2 S. 
The efficient synthesis on the surface and subse- 
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Fig. 13.— Calculated abundance, relative to hydrogen, of main sulfur species in ice mantle sublayers for 
CB 17 S and CB 26. The latter represents an significantly photoprocessed ice. 



































quent reactive desorption explain why H2S is the 
most abundant sulfur gas phase molecule for all 
modeled cases. 


Modeling results show that H2S and SO ices 
are produced during the freeze out epoch and con¬ 
tain most of the sulfur reservoir in interstellar ices. 
The reaction set for sulfur species is rather poor 
and does not include sulfurous and sulfuric acids 
(H2SO3 and H2SO4, respectively) that likely are 
very stable forms of sulfur in watery environment. 
Oxoacids and their derivatives are possible and 
non-det ectable forms of sulfur i n star-formation 
regions ( Kalvans fc Shmeldl 2010 ). Given the rich 
chemistry of sulfur and the radical nature of the 
SO molecule, it is unlikely that SO is a final, 
stable, and abundant form of S in interstellar or 
circumstellar ices. More likely, in terms of this 
model, it can be view ed as a representation of ox¬ 
idized sulfur (see also IScappini et aDl 2003 [l . 

Figure[T 51 shows an example of depth-dependent 
composition for sulfur species. H2S and SO are 
dominant at all depths, except for sublayer 3 at 
late times, where the abundance of NS is highly 
enhanced. Although this is probably because of 
the limited reaction network for sulfur, this illus¬ 
trates the oxidation of S during longer timescales. 
Oxidation products are SO, SO2, NS or perhaps, 
other species, not included in the network. Sulfur 
in the outer sublayer 1 is less susceptible to ox¬ 
idation thanks to a higher availability of atomic 
H. 


The high abundance of NS in sublayer 3 for 
CB 17 L, CB 17 S, CB 27 , and B 68 is an example 
how local conditions in the ice mantle may give 
rise to peculiar chemical features. NS is produced 
in all layers via the reactions N -|- S and NH -|- S. In 
sublayers 1 and 2 it is also efficiently destroyed via 
reactions with atomic C and N. Sublayer 3 is poor 
with CO and N2. Thus, C and N are unavailable, 
and NS accumulates after all remaining CO has 
been converted into CO2. 


tations. To improve this situation, we introduce 
two simple changes in the model - a modification 
to CO and H2CO hydrogenation activation ener¬ 
gies and, subsequently, a mild and short tempera¬ 
ture spike during the quiescent phase. The results 
of the three simulations are described in the fol¬ 
lowing sections. 

To provide an indication of the observed abun¬ 
dances of gas-phase COMs relative to the cal¬ 
culated values, the observed values have been 
marked in their corresponding figure panels. The 
observed values have no temporal or spatial con¬ 
straints with regard of the present model. The 
OD character of the simulation means that a qual¬ 
itative comparison is possible, only. We assume 
that the synthesis of an organic molecule is repro¬ 
duced with a degree of success, if its calculated ice 
abundance is not lower than its observed gas-phase 
abundance. 

5.1. COMs: Standard model 

5.1.1. The chemistry of COMs 

The ” Standard” model here is the one described 
in full in section [ 2 j Figures [ 14 ] and [ 15 ] show the 
calculated abundances in the center of the starless 
cores, along with known observational values. The 
species can be categorized in several overlapping 
groups: 

• those synthesized via the HCO radical, e.g., 
formic acid HCOOH, formaldehyde H2CO, 
methyl formate HCOOCH3, formamide 
NH2CHO, and acetaldehyde CH3CHO; 

• those synthesized via the CH3O radical, 
e.g., methanol; CH3OH, dimethly ether 
CH3OCH3, HCOOCH3, and as a subgroup, 
species whose synthesis involves CH3, which 
is a daughter of methanol. This subgroup 
includes CH3CHO, acetonitrile CH3CN, and 
methane CH4; 


5. Results: complex organic molecules 

When compared to observational evidence, the 
synthesis of organic species in the model can be 
characterized as moderately efficient. For some 
species, gas phase abundances are in temporal 
qualitative agreement with the observations, for 
others even ice abundances are well below expec- 


• carbon chain molecules, e.g., CH3CHO, 
propynone CH2CCO, ketene CH2CO, and 
acetylene C2H2; 

• species that are typically formed via atom 
addition on the surface, e.g., CH4, H2CO, 
HCO, HCOOH CH2CCO, CH2CO, C2H2, 
methyl imine CH2NH, cyanamide NH2CN, 
and CH3CN. A relatively high gas-phase 
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Fig. 14.— Calculated abundances, relative to hydrogen, of observed organic molecules, set No. 1. Gray 
lines: gas-phase abundance; black lines: abundance in ice (surface and three mantle sublayers). Contraction 
(Phase 1) ends and the stab l e core Pha se 2 begins at t?&es betw een 1.4 and 1.5Myr for all si mulations. Refer- 
ences: M - iMarcelino et aL ( 2005 ): O - Qberg et al.l ( 2010l) : B - iBacmann et al. ( 2012[ ): C - ICernicharo et al 
(l2012h : V - IVastel et all (l2014l) . 




























































































































































Fig. 15.— Calculated abundances, relative to H, of observed and other selected organic molecules, set No. 2. 
References as in figure [HI 
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Fig. 17. 
CB 17 L 





























abundance during the freeze-out epoch is 
characteristic for these species. 

The hydrogenation of CO and H 2 CO are the 
two most important steps that regulate the general 
production efficiency for many oxygen-containing 
COMs in bulk ice. Figures [16] and [17] show the rel¬ 
ative abundances of organic species in the four ice 
layers considered in the models for CB 17 L. This 
core can be regarded as a ‘median’ case among the 
five models, without the extremes in ice composi¬ 
tion shown by CB 17 S and B 68 on one hand, and 
CB 26 and CB 27 on the other hand lsection f3.ll) . 

Most of the species with abundances shown in 
figures fTCl and fTTl have their peak abundances high¬ 
est in the outermost sublayer 1. This indicates 
that even a partial ice sublimation might be suffi¬ 
cient to reproduce the observed gas-phase abun¬ 
dances of many COMs. In sublayer 1, CO is 
available in vast quantities at any integration time 
points, while the necessary H atoms are produced 
in H 2 O or H 2 photodissociation. With the present 
OD model it is not possible to determine the per¬ 
centage of ice mass that should be sublimated. If 
the modeled ices are assumed to be representative, 
this proportion would be ~ 1 — 10 %. 


5.1.2. Comparison with observations 


Figures [14] and [15] show that the model pre¬ 
dicts relatively high gas-phase abundance for some 
species in the HCO group - formic acid HCOOH, 
formaldehyde H 2 CO, and formamide NH 2 CHO. 
Although the observed gas-phase abundances are 
reached for a short period of time, this occurs dur¬ 
ing the infall Phase 1 and likely cannot be treated 
as an agreement with observations of stable star¬ 
less cores. 

For the CH 3 O group, the calculated abun¬ 
dances are significantly lower than those ob¬ 
served in interstellar conditions. Even assum¬ 
ing instant ice sublimation, observed gas-phase 
abundances cannot be reproduced for methyl 
formate HCOOCH 3 , dimethyl ether CH 3 OCH 3 , 
and methanol CH 3 OH. This result is similar 
to t hat of earlier prestella r core models with- 


out ( Garrod fc Herbslj 20061 : Garrod et al 


and with bulk ice photochemistry (iGarrod 


2008) 

20131 


Paper iJ ). It is this discrepancy that prompts us 
to proceed with a more detailed investigation on 
the chemistry of COMs. 


5.2. The Compromise model 

5.2.1. Changes relative to the Standard model 


In the present surface reaction network ([Garrod et al 


20081 ) the activation energy barriers Ea for reac¬ 


tions H-I-CO and H-I-H 2 CO are higher than 2000K, 
which means that they are inefficient for cold core 
conditions. Experimental evide nce points to much 


lower values for Ea, 300-600K ( Awad et al. 2005t 


Euchs et"^ 2 OO 9 I) . However, putting such low Ea 


values into the present model results in that al¬ 
most all carbon in ice is in the form of methanol, 
which is unrealistic. Higher Ea is especially re¬ 
quired for the CO-I-H reaction. 

We present results of model calculations that 
use the experimental, low Ea value for formalde¬ 
hyde hydrogenation. Eor CO hydrogenation Ea, 
an empirical compromise value was used, which 
lie s between the low experimentally detected Ea 
of Fuchs et ah ( 2009) and t he hig h (2500K) value 
proposed by Garrod et al.l ( 2008 ). Activation en¬ 
ergies for the three important CO hydrogenation 
reactions have been summarized in Table H We 
dub this the ” Compromise model” for further ref¬ 
erence. 


5.2.2. Comparison with observations 

Figure [TH] shows species, whose abundance is 
substantially affected with the changes in acti¬ 
vation energies for reactions 17, 18, and 19 in 
Table [H Comparison with the Standard model 
results (figures [14] and [15]) reveals that the ice 
abundance of species that are formed via CH 3 O 
is increased by up to six orders of magnitude 
(CH 3 OCH 3 in B 68 ). Importantly, the ice abun¬ 
dances of CH 3 OCH 3 , CH 3 OH, and CH 3 O are now 
above or similar to the detected gas phase values. 

The synthesis of methyl formate HCOOCH 3 
depends on both HCO and CH 3 O radicals. We 
found it impossible to achieve a sufficiently high 
ice abundance for HCOOCH 3 with changes in 
Ea for reactions 17, 18, and 19. Thus, the cal¬ 
culated relative ice abundance for this molecule 
remains approximately four orders of magni- 


tude below the obser 
are between 10 “^^ ai 

ved gas-phase values, which 

id 10 ® (Oberff et al. 

201(1 

Bacmann et al. 2012: 

Cernicharo et al. 2012 

)■ 


Methyl formate is also inefficiently produced in 
several other models considering either starless or 
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Fig. 18.— Calculated abundances, relative to H, of selected organic molecules for the Compromise model. 
Gray lines: gas-phase abundance; black lines: abundance in ice (surface and three mantle sublayers). For 
references see the caption of figure [TH 32 









































































































































Table 6 : Changes in reaction data for the compromise model. 


No. 

Reaction 


Standard® 

Activation barrier 
Fuchs et al. 12009) 

K 

Compromise 

17 

H-^CO ^ 

► HCO 

2500 

390-520 

1600 

18 

H-HH 2 CO 

CH 3 O 

2100 

415-490 

415 

19 

H-HH 2 CO 

CH 2 OH 

2500 

415-490 

415 


“The initial values used in this paper 


preste llar cores (lGa rrodll20 13t IVasvunin fc Herbst 
I2013al f9l. Chang fc Herbsd ( 2014 1 are able to effi¬ 
ciently synthesize HCOOCH 3 with a Monte-Carlo 
model, but their results also show a likely excess 
of CH 4 and a shortage of CO 2 - discrepancies gen¬ 
erally not observed in the abovementioned other 
models. These results may indicate that the in¬ 
efficient synthesis of HCOOCH 3 is not an arti¬ 
fact, caused by the limits of the present study. 
A higher temperature may be necessary for an ef- 
ficient synthesis of HC OOCH 3 , as suggested by 
ICarrod fc Herbstl ( 2006l l. 


5.2.3. Chemistry differences with the Standard 
model 

The most significant abundance increase for the 
CH 3 O group is observed in the B 68 model. First, 
the lower temperatures in this core mean a low 
mobility for atomic H in ice. Second, the abun¬ 
dant CO is now more easily hydrogenated, thanks 
to the lower Ea and CO soaks up much of the 
available atomic hydrogen in the mantle. These 
two factors result an accumulation of radicals so 
that they have greater chance to react with other 
multi-atom species. For the present discussion on 
COMs, the most important radicals are CH 3 O and 
HCO. The effect is especially visible for B 68 and 
CB 17 S. 

For species that are formed via HCO, the sit¬ 
uation is less clear because formaldehyde H 2 CO 
is consumed by reactions 18 and 19. In cores 
CB 17 L and CB 27, formaldehyde has its maxi¬ 
mum ice abundance decreased below the detected 
gas-phase values. A similar situation can be ob¬ 
served for formic acid HCOOH, directly related to 
HCO. 

Figure [19] shows that in the compromise model 
the abundances of COMs are increased in all sub¬ 


layers (cf. figures [T6l and flTll. The main part of 
their synthesis still occurs in the outer sublayer 1 . 
Species whose synthesis involve CH 3 O and CH 3 , 
do not reach a steady abundance until the end 
of the simulation, mainly because they can be re¬ 
garded as the end products of CO hydrogenation 
sequence. 

5.3. The synthesis of COMs in ice in a 
temporal -warm-up event 

The discussion in the above section 15.2.21 shows 
that for some COMs (notably, methyl formate) 
the calculated abundances are nowhere near their 
observed values. To test how a mild and temporal 
warm-up event can affect the ice abundances of 
COMs, we used the B 68 Compromise model that 
undergoes a temperature increase to 20K for Ikyr 
at t = 1.6Myr. 

5.3.1. Chemistry differences with previous mod¬ 
els 

Figures |20| and [21] show that the tempera¬ 
ture spike positively affects the production of a 
number of solid species - HCOOCH 3 , CH 2 CCO, 
CH 2 CO, NH 2 CHO, C 2 H 2 , CH 4 . Notably, the 
gas phase abundance is substantially increased for 
species synthesized on the surface - CH 2 CCO, 
H 2 CO, NH 2 CN, C 2 H 2 , CH 3 CN, CH 2 NH, and 
CH 4 , thanks to reactive desorption. 

For most species, the changes in ice abundance 
are within one order of magnitude. However, 
for methyl formate HCOOCH 3 and formamide 
NH 2 CHO an increase by approximately five and 
three orders of magnitude is observed, respec¬ 
tively. We therefore explore deeper the tempera¬ 
ture spike-induced chemistry of these two species. 
Initiated at 1.6Myr, the rise of abundances actu¬ 
ally occurs over a period of lOOkyr after the tem- 
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Fig. 20.— Comparison of calculated COM abundances, relative to H, for the B 68 core between the Standard 
model (st.), Compromise model (c.), and Compromise model with a warm temperature spike (c.w.), set 
No. 1. The two latter models begin to show differences only after 1.6Myr, when the temperature spike 
occurs. References as in figure [HI 
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Fig. 21.— Comparison of calculated COM abundances, relative to H, for the B 68 core between the standard 
model (st.), compromise model (c.), and compromise model with a warm temperature spike (c.w.) at 1.6Myr, 
set No. 2. References as in figure [TH 


35 








































































































Fig. 19.— Calculated abundances, relative to H, 
in ice sublayers for selected organic molecules, 
Compromise model, CB 17 L core. 


perature spike. Figure shows the abundances 
of NH 2 CHO and HCOOCH 3 in ice sublayers. 

Methyl formate production occurs mainly in 
the middle sublayer 2. This is unlike most other 
COMs in the Standard and Compromise models 
that are synthesized in the outer sublayer 1 (fig¬ 
ure [HI). Sublayer 2 is the only one that is rich 
in both, CO and H 2 O. The former is the source 
species for HCO and CH 3 O radicals that combine 
into HCOOCH 3 . Meanwhile, the photodissocia¬ 
tion of water is the principal source of atomic hy¬ 
drogen, necessary for the hydrogenation sequence 
CO ^ HCO ^ H 2 CO ^ CH30;CH20H. 
Thus, CO hydrogenation by the mobile H atoms 
at the 20K spike, followed by combination of HCO 
and CH 3 O Ri^lOOkyr after the spike, is the main 
mechanism for the formation of HCOOCH 3 . 

The case of formamide is different - NH 2 CHO 
reaches similar abundance in all three sublayers 
after the 20K spike. This is because in none of the 
sublayers its main parent species - CO and NH 3 - 
both have a high abundance. While NH 3 is mainly 
concentrated in sublayer 3, CO is abundant only 
in sublayers 1 and 2 . 

5.3.2. Discussion on temperature spike model re¬ 
sults 

The ice abundance increases occur because dif¬ 
fusion, proximity, and activation energy barriers 
are more effectively overcome at 20K. The sub¬ 
sequent reactions temporary reduce of ice abun¬ 
dances for the radical species HCO, CH 3 O, and 
CH 2 OH. The results shown in figures [HI and HT] 



Fig. 22.— Calculated abundances, relative to H, 
in ice sublayers for species most affected by a Ikyr 
20K temperature spike in the Compromise model 
for the B 68 core. 
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suggest that even such a small and short warm-up 
period is sufficient to qualitatively reproduce the 
abundance of methyl formate in ice and that the 
ice abundances of other organic species do not be¬ 
come disbalanced. The single-point OD approach 
used in the model prevents a quantitative compar¬ 
ison between calculation results and observations. 


We suggest two likely causes for a (temporal) 
temperature increase for a portion of ices in a star¬ 
less or prestellar core. First, the gas parcel can 
be exposed to the interstellar radiation because 



also result in mild heating of the gas. Both of these 
mechanisms also offer a hypothesis for the desorp¬ 
tion mechanism - either photodesorption or colli- 
sional sputtering of the weakly-bound, non-polar 
ice outer layers. These considerations are sup¬ 
ported by t he turbulent motions observed in star¬ 
less cores (iRedman et al.l 2006t Levshakov et al 
20141 : Steinacker et al.l 2014 1. 


6. Summary 


A model that describes ices in a detailed way 
was used to investigate ice chemistry in five dif¬ 
ferent examples of starless cores. While the im¬ 
portance of subsurface ice chemistry has been 
noted in previous studies (ICupoen fc Herbstll2007 : 
iKalvans fc Shmeld 2010 l : Icarrodl 2013[l . the most 
general conclusion from the present study is that 
chemical processes in ice depend strongly on the 
particular (sub)layer, where the molecules in con¬ 
sideration are located. This arises because of vari¬ 
ations in abundances of major ice components 
(H 2 O, CO, CO 2 , N 2 , NH 3 ) in different layers on 
the same grains. This conclusion is in line with 
the results of a prestellar core model ( Paper Ilh 
and a pseudo-time depende nt Monte Carlo model 
by ( Chang fc Herbst l2014h. Sing l e sub l ayer ap¬ 


proach (IKalvans fc Shmeldl 120101 l2013t iGarrod 


I 2 OI 3 I: iBelloche et al.ll20ll ~ is not suitable for rep¬ 
resenting subsurface chemistry of interstellar ices. 
This is especially so if the synthesis of minor 
species is studied. 


The H 2 O and CO ices in dark cores are pho- 
toprocessed to CO 2 via sequence O- This sets 
limits to core lifetimes, because current observa¬ 


tions indicate that the C 02 : H 2 0 abundance rati o 


in ice likely is 44% or less (jBoogert et al.N201lfl . 


The maximum lifetime for a quiescent and stable 
molecular core with constant physical conditions 
(Phase 2) inferred this way in the present model is 
650kyr (CB 17 S core). It is hard to pinpoint a pre¬ 
cise value because of the limitations in the physi¬ 
cal model, and uncertainties regarding the turbu¬ 
lence in dark cores. Taking this into account we 
can state that the model indicates dark core life- 
time s of <lMvr. The low er C 02 :H 2 0 ratio in Tau¬ 
rus (|Whittet et al.ll2007^ may indicate that these 
starless cores may be younger than those in most 
other mole cular clouds, observ ed bv iBoogert et al 
( 2011 1 and Oberg et al. ( 2011 1 . 


The H 2 0:C0 ice abundance ratio cannot be 
used as reliable indicator of core age because CO 
resides mostly on the outer surface of the grains. 
There, it can be rather easily lost to the gas phase 
or converted into other species in mild heating 
events. Such processes can be induced by expo¬ 
sure to interstellar photons or low-velocity shocks 
within the core, as discussed below. 


Because of more efficient desorption, the present 
model predicts shorter and higher abundance 
peaks during the cloud collapse phase (or ice for¬ 
mation epoch) for a number of species produced 
in dense gas or on the surface (section [331) • These 
include H 2 O, CO 2 , NH 3 , H 2 O 2 , H 2 S, SO, SO 2 , 
OCS, HCOOH, H 2 CO, CH 2 CO, CH 4 , O 2 , H 2 S 2 , 
N 2 etc. 


We find that O 2 and H 2 O 2 are the main oxygen 
reservoirs in ice, except water and carbon oxides. 
O 2 forms in the gas and on the surface, and is 
being converted into the H 2 O 2 via photprocess¬ 
ing in subsurface ice. It was suggested that H 2 O 2 
and O 2 H reach their gas-phase abundance peaks 
during the freeze out epoch, thanks to an active 
molecular exchange between the gas and the sur¬ 
face. The principal cause of desorption, necessary 
for this exchange, is irradiation by interstellar pho¬ 
tons. 


Calculations show that the HNCiHCN ratio in 
ice can be as high as 1.75 (in CB 27 core model), 
while it remains close to or below unity in the 
gas phase. Sulfur in ice is roughly equally divided 
between reduced (H 2 S) and oxidized (SO) forms. 
Subsurface processing of these sulfur species gives 
rise to OCS, NS, and H 2 CS, indicating an inter¬ 
esting diversity for the sulfur chemistry. 
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Investigating the chemistry of COMs, we found 
that lower activation energies - more consistent 
with experimentally detected values - for hydro¬ 
genation reactions of CO and H 2 CO may help ex¬ 
plaining the observations of COMs in the interstel¬ 
lar medium. A mild and temporal warm-up event 
helps to produce an abundance of methyl formate 
in ice that is higher and more consistent with gas- 
phase HCOOCH 3 observations. Such an event 
also substantially increases the ice abundance of 
formamide. 


Desorption by interstellar photons and ice sput¬ 
tering in grain collisions may transport the COMs 
to the gas phase. Both of these mechanisms re¬ 
quire either turbulence within the cores or external 
influence. These findings are probably supported 
by the fact that two of th e dense cores with de- 
tected CO Ms - L1689b (Bacmaiin et ahl 20121 
and Bl-b (Marcelino et al.l 20051: Oberg et al 


20101 : Cernicharo et al. 2012ll - have nearby star- 


form ation regions (IKirk et al.ll2007t lOber g et al. 
2 OIOII . while the third, LI544 ( Vastel et al.l 201 jb 


is kno wn to be turbulent and on the verge of co l- 
lapse ( Caselli et al. 2002 : Tatematsu et al. 2014h . 
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